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We study the scaling properties of heat transfer Nu in turbulent thermal convection at large 
Prandtl number Pr using a quasi-linear theory. We show that two regimes arise, depending on the 
Reynolds number Re. At low Reynolds number, NuPr~~ x ^ 2 and Re are a function of RaPr~ 3 ^ 2 . 
At large Reynolds number NuPr 1 ^ 3 and RePr are function only of RaPr 2 ^ 3 (within logarithmic 
corrections). In practice, since Nu is always close to i?a 1//3 , this corresponds to a much weaker 
dependence of the heat transfer in the Prandtl number at low Reynolds number than at large 
Reynolds number. This difference may solve an existing controversy between measurements in SF6 
(large Re) and in alcohol/water (lower Re). We link these regimes with a possible global bifurcation 
in the turbulent mean flow. We further show how a scaling theory could be used to describe these 
two regimes through a single universal function. This function presents a bimodal character for 
intermediate range of Reynolds number. We explain this bimodality in term of two dissipation 
regimes, one in which fluctuation dominate, and one in which mean flow dominates. Altogether, 
our results provide a six parameters fit of the curve Nu(Ra, Pr) which may be used to describe all 
measurements at Pr > 0.7. 

PACS numbers: 47.27 -i Turbulent flows, convection and heat transfer - 47.27.Eq Turbulence simulation and 
modeling - 47.27.Te Convection and heat transfer 



I. MOTIVATION AND OBJECTIVES 

In September 2001, a workshop on "high Rayleigh 
number convection" was held at Illmenau (Germany). 
During the course of this workshop, a controversy 
arose about the scaling and magnitude of the heat 
transport, in large Prandtl number experiments. Two 
groups, one in Santa Barbara pQ, and one in Hong 
Kong [21 [3] , reported measurements performed in wa- 
ter and organic fluids (alcohols) at Prandtl numbers 
ranging approximately Pr — 4 to Pr = 1300 and 
in cell of aspect ratio unity. In these experiments, 
the non-dimensional heat transfer between the bot- 
tom and top, Nu, appears to depend only weakly on 
the Prandtl number, and to increase with the Rayleigh 
number Ra, like: 



A third group of people, located in Rehovot, reported 
measurements performed with SF G gas, near critical 
point. By varying the distance to critical point, one 
can vary the Prandtl number between Pr — 1.4 up to 
Pr = 36, with relatively small non-Boussinesq effects 
up to Ra = 10 15 [I]. The heat transfer measurements 
in these experiments revealed a much stronger Prandtl 
number dependence, and appeared to be about 1.5 
larger in magnitude than in the water/alcohol exper- 
iments in the same range of Rayleigh numbers: 

Nu = 0.51Ra°- 29 Pr-° 15 . (2) 

Since the aspect ratio in the SFq experiments was 
comparable (0.72 v.s. 1) to the water/alcohol experi- 
ments, there was no clear explanation of the discrep- 
ancy between the two type of measurements. Com- 
pressibility effects due to the vicinity of the critical 



point were mentioned as a possibility, even though 
simple estimates seem to rule out their influence in 
the range of Rayleigh number explored [4] . 

Since then, the Rehovot group has improved its 
method of determination of the Nusselt number, 
which previously gave only precise estimate of rela- 
tive changes, not absolute ones. The improvement 
results in a better determination of the prefactor of 
the law |2]) which appears to be decreased by a factor 
2 with respect to the old measurements (Steinberg, 
private communication). This solves the discrepancy 
in magnitude between the Nusselt measurents in al- 
cohol/water and in SF6. But the discrepancy in the 
Prandtl dependence still remains. 

Clearly, in this type of controversy, it would be help- 
ful to rely on theoretical results to guide our intuition. 
The water/alcohol experiments appear to be in good 
agreement with a recent theory by Grossmann and 
Lohse |6] , which predicts a very weak decrease of 
the Nusselt number with the Prandtl number, due to a 
saturation of the viscous boundary layer (it increases 
with Pr). On the other hand, the SFq experiment is 
in excellent agreement with a theory by Shraiman and 
Siggia [7] , based on the existence of a turbulent bound- 
ary layer. Since both theories are mutually excluding, 
they can however not be used to solve the discrepancy. 
In a recent work [8HT0]. we have developed a quasi- 
linear theory of turbulent convection, which enables 
analytical predictions of the scaling laws of convection, 
including logarithmic corrections. The SFq experi- 
ment appeared to fit nicely into one scaling regime, 
in which the energy dissipation is dominated by the 
fluctuations, while the heat dissipation is dominated 
by the mean flow. The theory of Shraiman and Siggia 
also appeared as a special case of this theory, in which 
both the energy and the temperature dissipation are 
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dominated by mean quantities. So there are actually 
two regimes which could explain the SFq experiment. 
It is then natural to wonder first which regime is the 
more suitable for explaining the SFq and whether the 
quasi-linear theory can be used to propose solutions 
to the water/alcohol-ST 7 ^ controversy. This is the pur- 
pose of the present letter. 



II. THE QUASI-LINEAR THEORY IN 
SHORT: IDEAS, WEAKNESSES, 
CONSTRAINTS 



In this context, it is interesting to summarize the 
main ideas behind the quasi-linear theory. This theory 
is first used to compute analytically, from the Boussi- 
nesq equations, the mean and fluctuating velocity and 
heat profile within the turbulent boundary layer as a 
function of the Rayleigh and Prandtl numbers. Un- 
der the assumption that the energy and temperature 
dissipation are controlled by the turbulent boundary 
layer, these results can then be used to evaluate them 
analytically as a function of the Rayleigh and Prandtl 
number. Because this evaluation involves vertical in- 
tegration of the local dissipation over the height of the 
turbulent boundary layer, the final result also becomes 
a function of the turbulent boundary layer scale. In 
the quasi-linear approximation, this scale can be com- 
puted by assuming that the velocity or thermal fluctu- 
ations are passively advected by the mean flow [7J |S] • 
However, because this mean flow a priori originates 
from the bulk (it is external to the turbulent bound- 
ary layer), it cannot be analytically computed and 
has to be taken as an external parameter. The num- 
ber of degree of freedom can be greatly reduced if one 
assumes an algebraic profile U = z e for this mean 
flow. Indeed, by matching this mean flow with the 
analytic expression within the boundary layer, the pa- 
rameter e becomes the only free relevant parameter of 
the problem. In [10], we used velocity profiles given 
by numerical simulations by Kerr [TTJ fTJ] and Verz- 
icco and Camussi [T3] to postulate that for Pr < 0.7, 
the parameter e = 0, while for Pr > 0.7, its value 
is e = —1/2. This difference then provides a natural 
explanation of the difference between low and large 
Prandtl number measurements. Also, because bulk 
circulation can be severely affected by the geometry 
of the container (via the aspect ratio) , it introduces a 
non-trivial dependence in the geometry which is diffi- 
cult to quantify. In this regard, it is a weak point of 
the theory. 

On the other hand, we show below that the theory 
offers possibilities to connect different parameters of 
the convective cell independently of the bulk circula- 
tion, and depending only on the nature of the process 
which drives the dissipation (i.e. whether it is the 
mean flow or the fluctuations). Also, the hypothe- 
sis about the bulk flow contains several possibilities 



of direct check of the theory, via comparison between 
boundary layers length scales, and velocity profiles. 
Indeed, it can be shown [TU] that the exponent of the 
velocity profile and the length scale of the boundary 
layers are related through: 
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Here, u T = \J—vd z U x is the friction velocity, and 6(e) 
is equal to 1 if e = and zero elsewhere. All quanti- 
ties have been non-dimensionalized by the vertical gap 
and the thermal diffusivity. The subscript T stands 
for thermal, and V stands for velocity boundary layer. 
The dependence of e in Xblt and Xblv has been 
omitted for simpler notations. Also, the logarithmic 
and constant velocity regime have been lumped into 
the single notation U(X) <~ u T /(l/lnA + 1//) which 
patches the two consecutive characteristic behaviors. 
Note that in these units, u T /Pr = Re*, the Reynolds 
number based on the friction velocity |19j. Equation 
(|3|) then encompasses well-known relations. For e = 0, 

— —1/2 

one finds Xblv ~ Re* > like in a typical (logarith- 
mic, i.e. e = 0) boundary layer. For e = 1, like in a 
laminar layer, one finds Xblv ~ Re 
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This is the 



scaling originally proposed by Shraiman and Siggia 

ra- 



in. APPLICATION TO LARGE PRANDTL 
NUMBER CONVECTION 

A. Reynolds number 

What would be the situation in large Prandtl num- 
ber convection? Because v ^> 1, it can be ex- 
pected that small-scale temperature fluctuations, if 
any, are mainly affected by the laminar part of the 
velocity field. Applying ([3]) with e = 1, one then 
expect that the thermal boundary layer scales like 

Xblt ^ Pr 1 / 3 ?v 2/ ' 3 ■ If the thermal fluctuations are 
so weak that only the mean temperature dominates 
the dissipation, then Xblt ~ 1/Nu [ID], fixing the 
link between Re* and Nu as: 



Re* 



Nu 3/2 Pr~ 1/2 . 



(4) 



These laws are quite general, and independent of the 
shape of the velocity in the bulk region. We thus 
expect them to be universal laws of the large Prandtl 
number regime. So, is it satisfied by the available 
measurements in both water/alcohol and SFq? 

In the SF6 case, there are no direct measurements 
available of the friction velocity Reynolds number 
Re*. Instead, Ashkenazi and Steinberg [1] provide 
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measurements of a Reynolds number based on the 
peak frequency f p from the velocity power spectra 
Re — Af p L 2 /v. It can be fitted by a power-law 
form Re = 2.6Ra 0A3±0 ' 02 Pr-°' 74±0 02 . Comparing 
the power-law form of Nu and Re, one then sees that 
the law Q is very well satisfied provided Re ~ Re*, 
as already stressed in [3]. 

In the alcohol case, the Hong-Kong group measured 
Nu = Ra ' 278 Pr~ 028 . Other direct measurements 
concern not the friction Reynolds number, but the 
large scale Reynolds number Re. Its variation with 
Pr and Ra has been summarized by the Hong Kong 
group as: Re = 1.09Ra OA3 Pr~ ' 76 [3 . Interestingly 
enough, this variation is very similar to the variation 
deduced in SF6 (with a smaller prefactor, resulting in 
smaller Reynolds numbers). However, given the very 
weak dependence of Nu with Pr, this measurements is 
clearly incompatible with Q , which shows that in the 
alcohol case, either the equality Re ~ Re* does not 
hold, or/and the general law Q is violated. In that 
respect, it is interesting to note that a direct mea- 
surements of Re* performed by the same group, in a 
water experiment at Pr = 7 for Ra between 10 7 to 
10 11 [14J provides additional support of the violations 
of this laws. Indeed, they found Re ~ Ra ' 5 , while 
Re* ~ Ra 33 , closer to Re* ~ Nu ~ Re 2 ! 3 than to 
Q. These interesting new scalings are reminiscent of 
measurements performed on the r.m.s. vertical veloc- 
ity fluctuations by the Rehovot group, which indicate 
that the corresponding Reynolds number obeys a scal- 
ing transition from Ra 0A3±0 " 2 at Pr = 27, 45 and 93, 
to a scaling Ra 34 at Pr = 190 [15]. 

The alcohol measurements show that our initial 
naive expectation gathered in (|4) has to be refined. 
For this, let us come back to the main hypothesis 
governing Q. It is based on the idea that the ther- 
mal boundary layer is determined by the balance be- 
tween the large scale advection of temperature fluc- 
tuations Ud9 and the heat dissipation d 2 9 ~ 0Nu~ 2 . 
Within the boundary layer, the flow is laminar, ori- 
ented in the longitudinal direction and obeys the fric- 
tion law U = PrRe 2 z, where z is the vertical co- 
ordinate [10] . To obtain a "typical" velocity advec- 
tion, one then sets z — Xblt = Pr~ x l 3 Re*~ 2 ^ to get 
Re 2 *Pr = Nu 3 . This argument does not take into ac- 
count the constraint that the typical size of the bound- 
ary layer cannot increase beyond the size of the box 
[6] . This effect is obtained when the whole cell is lam- 
inar, i.e. as the Reynolds number decreases. To take 
this into account, we may then change our estimate 
into U = PrRel/(A + Pr^ 3 Rel /3 ), resulting into the 
more general law: 
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experiment of Ashkenazi and Steinberg) , one observes 
a transition from Re* ~ Nu 3 l 2 Pr~ x l 2 (i.e. the naive 
law Q) to Re* ~ NuPr~ x l 2 . 

Large scale circulations are easier to measure than 
friction velocities. Its is therefore interesting to con- 
nect Re and Re* to be able to exploit and verify ([5|. 
For this, we match the large scale circulation velocity 
RePr to the velocity vbl at the top of the boundary 
layer. At large Reynolds number, the boundary layer 
is turbulent and Vbl = Re*Pr = RePr (within log- 
arithmic corrections). At smaller Reynolds number, 
the boundary layer is laminar and vbl — Re 2 Xblt- 
The link between Re and Re* then depends on the 
shape of the velocity in the bulk flow (through the 
value of e). However, we may note that this length- 
scale is also subject to the condition that it cannot 
increase beyond the size of the box [6] (a situation 
occurring at very low Reynolds numbers). This re- 
marks allows us to build a general law resulting from 
a patching between the large Reynolds number case 
and the very low Reynolds number case as: 
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1 + Re*/Re c 



(6) 



where Cq is a constant depending a priori on the aspect 
ratio. Decreasing Re* (e.g. by increasing Pr as in the 



where c\ is a non-universal constant and Re c is a crit- 
ical Reynolds number. A rough estimate of this two 
parameters can be made using the alcohol data of Xia 
and collaborators. The available measurements un- 
fortunately only concern Nu and Re as a function 
of Ra and Pr. However, in this low Reynolds num- 
ber case, we expect R* ~ NuPr~ x / 2 . Using this, 
we may then find c\ and Re c by a two parameter fit 
of Re versus NuPr^ 1 / 2 . Fig. [I] shows a compari- 
son of the measured Re and the fit using ([6]) with 
C\ = 7, Re c = 17. The agreement is quite satis- 
factory. A further test of this fit can be made by 
comparing the effective boundary layer lentgh scale 
Xblv = Re/Rel = 1/(1 + Re*/Re c ) with the direct 
measurement of Xian and Xia [3HH] in water (Pr = 7) 
and alcohol, giving Xblv ~ Ra~ 016 Pr ' 24 . One sees 
in Fig. [2] that the agreement is satisfactory. Finally, 
one can invert ^ to obtain Re* as a function of Re 
in the SF6 experiment, using the values of the pa- 
rameters calibrated from the alcohol experiment. The 
results is shown in Fig. [3] It indeed give Re ~ Re*, 
as expected in this large Reynolds number case. 

The combination of ^ and |5]) predicts that, for 
a given (large) Pr, large Reynolds number leads to 
Re* — Re — Nu 3 l 2 Pr~ Y l 2 , while at lower Reynolds 
number, Re = Re* = Nu 2 / Pr. In between, there is a 
transitional regime in which Re varies smoothly with 
Re*. For the range of Reynolds numbers explored 
by the Hong-Kong group, we checked that this corre- 
sponds to an apparent scaling behavior Re = Re 2 J 3 = 
Nu 3 / 2 Pr~ 3 / 4 . The first case is representative of the 
data by Ashkenazi and Steinberg, while the second 
is representative of the data by Xia and collabora- 
tors, and Ahlers and collaborators. It seems therefore 
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FIG. 1: Reynolds number versus Nusselt number in the 
alcohol experiment [2]). The symbols are the experimental 
measurements. The line is the fit predicted by the model 
(eq. Q with ci = 7 and Re c = 17). 
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FIG. 3: Friction Reynolds number in the SF6 experiment. 
The symbols is the estimated friction Reynolds number 
derived using the theoretical formula (eq. (|6| with ci = 7 
and Re c = 17). The line is a fit Re* ~ Re, a behavior- 
predicted by the theory. 
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certain range of Reynolds number, the mean compo- 
nent of the turbulent flow undergoes a global bifurca- 
tion with three distinct flow configurations, character- 
ized by different transport properties, and subject to a 
strong hysteresis [16] . Is one of these two explanations 
sufficient to explain the difference in Prandtl number 
dependence observed? To answer this question, we 
now focused on the second ingredient of [SJ [71 [TO] : the 
energy dissipation, enabling to link Re, and Ra and 
Nu. 



B. Dissipation 



FIG. 2: Boundary layer scale in the alcohol experiment 
at 3 < Pr < 10). The symbols are the model prediction 
5 V = 1/{1 + Re/Yl). The line is the fit S v = Ra~° 16 'Pr - 24 
proposed by [2]. 



that we have identified a possible cause of discrep- 
ancy between the two groups: the difference in mag- 
nitude in the Reynolds number, resulting in a different 
regime explored. This difference of Reynolds number, 
observed in experiments with similar aspect ratios, 
Prandtl numbers and Rayleigh numbers is surprising 
at first sight. A first explanation is that compressibil- 
ity effects do play a role. A second explanation is that 
for this range of parameters, there are two possible 
solutions, with distinct transport properties (leading 
to different large-scale circulations and different Nus- 
selt numbers). Note that a similar phenomenon has 
recently be observed in the flow between two coax- 
ial rotating cylinders (the von Karman flow). In a 



In [TO], we have shown that the contribution to 
the energy dissipation scales like u^Xblv I Pr 
Pr 3 Re^\BLV for a laminar mean flow, 



,3 _ 



Pr 3 Re 



for a turbulent mean flow (larger Reynolds number), 
and like {PrNuf/ 2 Ra 1 / 2 ln\ B LV I '^blv for fluctua- 
tions. The first two contributions can be patched to- 
gether as a function of the Reynolds number in a way 
similar to what we did for the Reynolds number. Sum- 
ming the contributions of the mean flow and of the 
fluctuations (in a way similar to a procedure adopted 
by Grossmann and Lohse E] ) and recalling that this 
energy dissipation equals PrRa(Nu-l) [7J, we obtain 
in this framework a general relation linking Re* and 
the other parameters: 



PrRa{Nu-l) = c x Pr 



1 Z T*m +C2(PrNu)^Ra^ln\ BLV /X 2 BLV 
1 + Re*/Re c 



(7) 

where C\ and C2 are (a priori non universal) constants. 
Using (f5j, this relation can also be written in a more 



5 



illuminating form as: 



Ra 

p r 3/2 



Cl 



Nu 
Pr 1 / 2 



Re: 



1 + co 



Pr^ReJ 3 Nu 1 ' 2 Ra 1 ' 2 
+c 2 



1 + i?e*/J?e c 



p r l/4 p r 3/4 

(8) 

This relation can be used to determine interesting- 
scaling properties independent of the bulk flows or of 
the geometry. Indeed, we note that Xblv is a function 
of i?e* (see Further, we see from ^ that at low 
Reynolds number i?e* is a function only of NuPr^ 1 ! 2 . 
Considering ^ at low Reynolds, this means that in 
this regime NuPr^ 1 / 2 and Re are function only of 
RaPr~ 3 / 2 . This conclusion is independent of the fact 
that the dissipation is dominated by the mean flow or 
by the fluctuations. In practice, since Nu is close to 
Ra 1 / 3 and Re close to Ra 1 / 2 , this corresponds to a 
very weak dependence of Nu with Pr, and a depen- 
dence p r -0 ' 75 for Re. 

At large Reynolds number, i?e* is a function of 
NuPr" 1 / 3 and Xblv = Rel^ 3 [TU]. Inserting this into 
([8]), we see that at large Reynolds number, Re*Pr 
is a function only of RaPr~ x / 3 (within logarithmic 
corrections), or equivalently NuPr 1 / 3 is a function 
only of RaPr 2 / 3 . For Nu and Re close to a Ra 1 / 3 
and Ra 1 / 2 dependence, this means a decrease of Nu 
like Pr -1 / 9 (within logarithmic corrections) and of 
Re like Pr~ 2 / 3 . The decrease in Nu in that regime is 
therefore stronger than in the low Reynolds number 
regime, and qualitatively explains the difference be- 
tween the experiments in SF6 and in alcohol/ water. 
For a quantitative agreement, one would need to eval- 
uate the constants ci,..C2 appearing in ([8| and ([5]). 
One possibility would be to perform a non-linear fit 
of the experimental data points, like [BJ. In the sequel, 
we would like to explore another possibility, based on 
a collapse method inspired from finite size scaling in 
statistical mechanics. 



C. A scaling approach 

Equation Q, valid in both the large and low 
Reynolds number regime, provides an interesting im- 
plicit equation for the variable y = \J Ra/(NuPr): 



that as P — ¥ oo, i?e* 
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a 2 1 + cpPr^ReV 3 , ; x 
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y * 1 + Re*/Re c 
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(9) 

The function y{R Sf ,Pr) can be found once the func- 
tional shape for Xblv is provided. In the sequel, we 
shall adopt 



^blv = 1 + C 3-R* 



1/3 



(10) 



patching between the large Reynolds number case, 
and the saturation expected at very low Reynolds 
number. 

Our goal is now to determine the constants C\ , c 2 , C3 
and Re c from the available data. For this, we remark 



Note that given Q, this corresponds to Nu ~ Ra 1 / 3 , 
a limiting law also found by (0). In addition, exper- 
imental data indicate that in the range of parameters 
usually explored, the behavior of Nu with Ra and Pr 
follows approximate power laws. This suggests that it 
should be possible to approximate the behavior of y 
at large but finite Prandtl number via a scaling shape 
(reminiscent of finite size scaling in statistical mechan- 
ics): 



y 
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(12) 



where a is a scaling exponent (equivalent to a critical 
exponent) and F is a (universal) scaling function. In 
this approximation, the two parameters (a and F) can 
be found experimentally from a collapse procedure on 
the data taken at different Pr and i?e*. This pro- 
cedure is described below. We stress that since the 
scaling form is not exactly satisfied by the theoreti- 
cal formula Q, it will only provide an approximate 
experimental verification of ([9]). However, since this 
method uses data in different experimental set ups 
(different Pr and Re*), it is less liable to systematic 
experimental errors and thus, provides a more robust 
estimate of the constants c\, ...C3 than via a direct fit 
of the data. In addition, it provides a strong indepen- 
dent test of the theoretical formula ^ since it does 
not use its functional shape a priori. 



D. A collapse method 

In recent convection experiments, the quantity i?e* 
is usually not directly measured and we have to use 
only an estimate of it, based on Q. To minimize the 
error introduced by this procedure, we may rewrite 
(TT1) as: 



yPr a = G(Re*Pr a ) 



(13) 



with G(x) = yj cix 2 F{x). The determination of a and 
G is then made by trial and error, plotting yPr a as a 
function of Re*Pr a , and choosing a of that all curves, 
corresponding to different measurements at different 
Pr and Re, collapse onto a single curve. The resulting 
curve is the scaling function G. Fig. [IJpresent the best 
collapse obtained with a = 0.7 for the data collected 
by the groups of Ahlers (Santa Barbara), Xia (Hong- 
Kong) and Steinberg (Rehovot) [5D]. In the first two 
cases, we used i?e* — ^Nu 2 /Pr. In the third case, 
we used i?e* = 0. 17 yj Nu 3 /Pr. Interestingly enough, 
the scaling function seems to possess two branches 
at intermediate Reynolds number (one extending the 
Ahlers data, the other one extending the Steinberg 



G 



data). Note that among the data of Steinberg and his 
group, part of the data do fall onto the lowest branch 
corresponding to the Ahlers data. To check whether 
this is an artefact of our method, or whether it cor- 
responds to a real physical process (like bimodality of 
the Nusselt numbers, discussed in [17l [18], we ran an 
additional check. We present in Figure [5] the collapse 
obtained when the data of Chavanne et al [T7] have 
been include d, with t he Reynolds number estimated 
as Re* = \f 'Nu 2 / 'Pr. As can be seen, the univer- 
sal scaling function seems also to describe these data. 
Moreover, the range where our bimodality was de- 
tected appears to coincide with the bimodality range 
of Chavanne et al. Such a bi-modality could arise 
for example from a global bifurcation of the turbulent 
mean flow, like in the von Karman experiment (see 
Section III. A.) 
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FIG. 4: The universal scaling function. Best collapse of 
the data in alcohol [2], in water [l] and in SF6 [4], obtained 
by plotting V RaPr 0A Nu- 1 as a function of Re*Pr°.7. 
Note the bimodal regime appearing in intermediate range 
of Re t Pr°- 7 . 
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FIG. 5: Test of the bimodality. Same as Fig. 4, with 
the data of [T7] plotted as open symbols. The figure has 
been enlarged around the intermediate range of Re,Pr 07 
for more clarity. Notice the persistence of the bimodality, 
although it is slightly less apparent. 



The "lower branch", on the other hand, is well fitted 
by a "fluctuation contribution" : 



Vfl 
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(15) 

The two corresponding regimes (fluctuation or mean 
flow dominated) seem to be mutually excluding, ex- 
plaining the bimodality. We do not have an explana- 
tion for this. The quality of the fit can be checked a 
posteriori by solving simultaneously ( |15[ ) or ( 14 ) and 
Q for a given Ra and Pr, to find Nu(Ra, Pr)7 This 
consistency check is shown in Fig. [7]and|8]for the data 
of Ahlers, Xia and Steinberg. Overall, the quality of 
the fit is good. We expect improvement to be gained 
once direct measurements of Re* will be available. 



IV. SUMMARY 



E. Fits and consistency check 

This bimodality can be reproduced within the the- 
oretical model developed in Section III.B. In this 
model, we showed that the dissipation is the sum of 
two contribution: one from the mean flow (term pro- 
portional to ci), and one from the fluctuations (term 
proportional to C2). We show on Figure [6] that the 
upper branch of the universal function is well fitted 
by a "mean flow contribution": 
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We have shown that a discrepancy between heat 
flux in alcohol/water and in SF6, at large Prandtl 
numbers, could be explained by a Reynolds number 
effect and the existence of two different regimes, in- 
duced by a global bifurcation in the mean flow. We 
further showed how a scaling theory could be used to 
describe these two regimes through a single universal 
function. This function presents a bimodal charac- 
ter for intermediate range of Reynolds number. We 
explained this bimodality in term of two dissipation 
regimes, one in which fluctuation dominate, and one 
in which mean flow dominates. Altogether, our results 
provide a six parameters fit of the curve Nu(Ra, Pr) 
which may be use to describe all measurements at 
Pr > 0.7. This is one more parameter than the the- 
ory recently developed by Grossmann and Lohse [B], 
but it also fits the data by Steinberg and it may lead 
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to more precise estimates, as the fit is done simulta- 
neously over different experiments. 



nr 



FIG. 6: Fit of the universal function. The upper branch 
is fitted with formula ( 14 1 (solid line) , corresponding to 
a mean flow dominated dissipation. The "lower branch" 
is well fitted by a "fluctuation contribution" formula (151, 
(dotted line). 




FIG. 7: Nusselt versus Prandtl in the alcohol experiment. 
The open symbols are NuRa~ 25 for Ra = 10 9 ' 25 (circles) 
and Ra = 10 7 ' 25 (squares). The line are the theoretical 
predictions, using the fit ( |14[ ). 
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